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ABSTRACT 

The detection of high contrast companions at small angular separation appears feasible in conven- 
tional direct images using the self-calibration properties of interferometric observable quantities. The 
friendly notion of closure-phase, which is key to the recent observational successes of non-redundant 
aperture masking interferometry used with Adaptive Optics, appears to be one example of a wide 
family of observable quantities that are not contaminated by phase-noise. In the high-Strehl regime, 
soon to be available thanks to the coming generation of extreme Adaptive Optics systems on ground 
based telescopes, and already available from space, closure-phase like information can be extracted 
from any direct image, even taken with a redundant aperture. These new phase-noise immune ob- 
servable quantities, called kernel-phases, are determined a-priori from the knowledge of the geometry 
of the pupil only. Re-analysis of archive data acquired with the Hubble Space Telescope NICMOS 
instrument, using this new kernel-phase algorithm demonstrates the power of the method as it clearly 
detects and locates with milli-arcsecond precision a known companion to a star at angular separation 
less than the diffraction limit. 

Subject headings: techniques: high angular resolution, image processing, interferometric; methods: 
data analysis; stars: low-mass, close binaries 



1. PHASE IN THE FOURIER PLANE 

Only two parameters essentially determine whether a 
source is detectable during an observation: its bright- 
ness at the wavelength A of interest and the angular res- 
olution necessary to separate the source or feature from 
its direct environment. The angular resolution is ulti- 
mately constrained by the diffraction of the telescope, 
and astronomers usually follow the rule of thumb known 
as the Rayleigh criterion, stating that to be resolved, two 
sources need to be separated by 1.22A/D, where D is the 
diameter of the telescope used, to design their observa- 
tions. 

The development of optical interferometry has how- 
ever made this criterion obsolete: thanks to the exquisite 
level of calibration it permits, interferometry indeed 
makes it possible to detect sources or constraint the 
extent of features around objects at separations sig- 
nificantly smaller than the diffraction limit. Even at 
the scale of one single telescope, the results obtained 
with the technique known as non-redund ant masking 
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Adaptive Optics (AO) systems jlTuthill et alj|2005 [2006t 
Llovd et alj|2006t llreland et al.l T2008: Kr aus et al.ll200a 

Martinache et al1l2009l ) demonstrate the relevance of this 
technique for the detection of structures at small angu- 
lar separation, tha t would not be accessib le from conven- 
tional AO images (jRaiagopal et al.ll2004[) . 

Even if it only uses one single telescope, in order to 
reach this level of resolution, one however needs to ac- 
cept that the familiar product called "image" may not 
necessarily constitute the best final data-product. In- 
stead, when interested in high-angular resolution prop- 
erties of partially resolved objects, it is convenient to de- 
rive information not from the image itself, but from its 
Fourier-transform counterpart. This information, known 
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(NRM) interferometry, first , seeing- limited (|Haniff et aLl 
Il987t iRcadhca d et alJ ll988) and more recently used with 



as complex visibility, is extracted from the Fourier-plane, 
calibrated and then tested against a model of the ob- 
served object. 

In optical interferometry, this approach is often 
mandatory: the paucity of apertures (N ^2 — 5) and 
baselines make the content of a direct (Fizeau) image 
of limited value. Information rich images can be recon- 
structed after extraction of the complex visibility func- 
tion from the u, v plane, but only with a large (N > 10) 
number of apertures like in radio-interferometry, or after 
using image synthesis. The optical image reconstruction 
known as pupil densifi cation that is used in hypertele- 
scopes ()Labevrielll996f ) does provide an alternative, but 
again, only be comes compelling if a large number of aper- 
tures is used (|Labevrie et al.l 120081 ). But even when an 
image can be reconstructed from optical interferometry 
measurements , e.g. t he images of the binary Capella by 
iBaldwin et al.l (119961). the intens ity map of the surface 
of Altair bv iMonnier et all ()2007f ) or the spectacular im- 
ages of the disk eclipsing e Aurigae (Kl oppenborg et all 
l2010f ) , quantitative characteristics of the sources can only 
be deduced from the fit of the interferometric data by 
parametric models. In the case of a marginally resolved 
binary star, precise measurements of angular separation, 
orientation and contrast, with confidence intervals, de- 
duced from a model-fit of complex visibilities carry much 
more scientific value than an image of "blurry blobs" . 

Visibilities in the Fourier-plane are complex numbers, 
whose amplitude and phase are usually considered sepa- 
rately. This paper focuses on the treatment of the phase 
and ignores the amplitude. In general, the power con- 
tained at given spatial frequency is the result of the co- 
herent sum of R random phasors, with R a scalar coding 
the redundancy of the spatial frequency In the presence 
of residual optical path differences (OPD), this coherent 
sum of R random phasors loses the phase information 
and results in the formation of speckles in seeing-limited 
images with a visible/IR telescope. NRM- interferometry 
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solves this problem, by discarding light with a pupil mask 
designed so that each baseline is unique (i? = 1), which 
makes the extraction of the phase possible. 

The phases alone, being corrupted by residual OPDs, 
are of restricted interest. It is however possible to 
combine them t o form what is known as closure-phase 
(|Jennisonlll958D . that is the sum of three phases mea- 
sured by baselines forming a closed triangle. This re- 
markable interferome tric quan t ity (cf . the introduction 
to closure phase by iMonnierl ()2000D ) exhibits a com- 
pelling property: it rejects all residual pupil-plane phase 
errors. Moreover, because it is determined from the anal- 
ysis of the final science detector and not on a separate 
arm wavefront sensor, it is also immune to non-common 
path errors between the wavefront sensor and the science 
camera. 

Once extracted and calibrated, the closure-phases can 
then be compared to a parametric model, for instance 
of a binary star, to confirm or infirm the presence of 
a companion around a given source, while uncertainties 
provide contrast detection (i.e. sensit ivity) limits . This 
approach was successfully used by J Llovd et all 120061; 
Martinache et all l2007t llreland et a l.N2008l: iK raus et al.1 
2008HMartinache et al.ll2009t ). who typically report sen- 
sitivity of 5-6 magnitudes in the near infrared at separa- 
tions ranging from 0.5 to 4 X/D. 

This paper aims at generalizing the notion of closure- 
phase, and shows that closure-phase like quantities, i.e. 
sharing the same property of independence to pupil-plane 
phase errors, can be constructed even in the case of re- 
dundant apertures. 

2. KERNEL-PHASE 

2.1. Linear model 

Whether contiguous (i.e. single-dish) or not (i.e. in- 
terferometric), the 2D pupil of an imaging system can 
be discretized into a finite collection of N elementary 
sub-apertures. One of these elementary sub-apertures 
taken as zero-phase reference, the pupil-plane phase of a 
coherent point-like light source can be written as a N — 1- 
component vector ip. Given that the image, or interfer- 
ogram, of this source is sufficiently sampled, then in the 
Fourier plane (a.k.a. (u, w)-plane in interferometry) one 
will be able to sample up to M phases, where M is a func- 
tion of the pupil geometry only. For a non-redundant ar- 
ray made of N elementary sub-apertures, the number of 
sampled (u,v) phases is maximum M = (2)- The same 
number of sub-apertures organized in a redundant ar- 
ray, for instance following a regular grid, produces signif- 
icantly less distinct (u, v) sample points as each point re- 
ceives the contribution of several pairs of sub-apertures. 

In most cases, since each point receives the sum of sev- 
eral random phasors, both phase and amplitude are lost 
and cannot be simply retrieved: this results in the forma- 
tion of speckles. However, if the Strehl is high enough, 
the complex amplitude associated to the instrumental 
phase in one point of the pupil, tpk , can be approximated 
by e lipk »1 + iifk- Direct application of the approach is 
therefore for now, restricted to space-borne diffraction- 
limited optical and mid-IR telescopes like HST (cf. Sec- 
tion [3]), but should also prove relevant to the upcoming 
generation of extreme AO systems. 

Given that the proposed approximation holds, while 
observing a point source, the unknown (instrumental) 
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Fig. 1. — Iterative process for the determination of the transfer 
matrix A. The top row shows the sub-aperture of the full 2D pupil 
(circular aperture with 30% central obscuration) where a phase off- 
set is applied (three cases are represented) . The bottom row shows 
the resulting distribution of phase in the Fourier plane. The over- 
laid dashed-line circle in the bottom row marks the cutoff spatial 
frequency of the transfer function. 

phase distribution in the pupil ip can be related to the 
phases $ measured in the Fourier plane with a single 
linear operator. To build an intuitive understanding of 
this relation, let us consider the following scenarios: 

• If the phase is constant across the entire pupil, then 
none of the baselines formed by any pair of elemen- 
tary sub-apertures does record a phase difference, 
and the phase in the Fourier plane is zero every- 
where. 

• If a phase offset 8^ is added to one single sub- 
aperture, then each baseline involving this aperture 
records a phase difference, which is exactly ±5$. 
Figure [1] represents several such scenarios. 

• If the pupil-plane phase vector ip is completely ran- 
dom, each of the M samples in the Fourier-plane 
is then the average of R phase differences on the 
pupil, where R is the redundancy of the considered 
baseline. 

To reproduce this behavior, the following linear model 
will be used: 



$ = R 



A-ip, 



(1) 



where $ represents the M-component Fourier plane 
phase vector, R a M x M diagonal matrix whose di- 
agonal elements code the redundancy of the baselines, 
and A represents a M x (N — 1) transfer matrix, whose 
properties form the core of the discussion of this work. 
To be complete, the model should also include the phase 
information intrinsic to the observed source, represented 
by the term <&o that simply adds on top of the instru- 
mental phase. One can then multiply both sides of the 
equation by the matrix R so that it becomes: 

R$ = A-p + K-^o- (2) 

While R and A could have been merged into one sin- 
gle operator, there are intentionally kept distinct. The 
rationale for this choice so is that the left hand side of 
Eq. [21 i.e. the measurements, can be acquired by read- 
ing directly the imaginary part of the complex visibility. 
Given that the next (quadratic) term in the Taylor ex- 
pansion of e lv being real, this makes the approximation 
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valid to the third order in phase. This also makes A of 
striking aspect as it is then exclusively filled with values 
0, 1 or -1. 

If the matrix A were invertible, then the analysis of 
one unique focal plane image of a single star (case cor- 
responding to Eq. [IJ would be sufficient to determine 
the instrumental phase <p as seen from the detector, and 
drive an AO system and/or delay lines. Except for the 
special case of a non-redundant aperture, the problem 
is however known to be degenerate, despite the larger 
number of measures than unknowns (M > N — l). 2 

As demonstrated by the successes of NRM- 
interferometry, a complete characterization of the 
wavefront is not essential if one can determine observ- 
able quantities that are pupil-phase independent. The 
closure relations used in interferometry can be related 
to the operator A: these relations are simply linear 
combinations (modelized by an operator K) of rows of 
A that produce 0, forming something known as the left 
null space of A: 

K A = 0. (3) 

For a non-redundant array, each closure relation will 
fill a row of K with mostly zeroes, except in three po- 
sitions corresponding to the baselines forming a closing 
triangle, that will contain 1 or -1. These relations are 
however not the only possible ones, and less trivial com- 
binations, involving more than three rows at a time, can 
be constructed. The total number of independent rela- 
tions however remains unchanged and is only imposed by 
the geometry of the array. 

Although not impossible, finding the operator K "by 
hand" (i.e. finding a basis for the left null-space of A) 
for a redundant aperture is a tedious task, as the matrix 
A can get quite large. A very efficient way to do this 
is to calcutate the singular value decomposition (SVD) 
of A T . The SVD algorithm (jPress et all 120021 ). allows 
to decompose the now (N — 1) x M matrix A T as the 
product of a (N — i) x M column-orthogonal matrix U, a 
MxM diagonal matrix W with positive or zero elements 
(the so-called singular values) and the transpose of an 
M x M orthogonal matrix V: 

A T = U • W • V T . (4) 

One relevant property of the SVD is that it explicitly 
constructs orthonormal bases for both the null-space and 
the range of the matrix A T . Of particular interest here, 
are the columns of V that correspond to singular values 
equal to zero: these vectors form an orthonormal base 
for the null-space, also refered to as Kernel of A T , that 
is exactly what is needed to fill in the rows of K. 

If the observed target is not perfectly symmetric, and 
exhibits actual phase in formation (i.e. $o 0, see for 
instance [Monnieri (|2OOO0 ). Eq. [2]is required. Multiply- 
ing it with the left side operator K leads to a new se- 
ries of new phase-like quantities that are not contam- 
inated by instru mental phase, genera lizing the notion 
of closure-phase (|Baldwin et al.lll986|) on which N RM- 
interferome try from the ground (|Tuthill et al.ll2000T) and 
from space ( Sivaramakrishnan et al. 2009) rely entirely. 




Fig. 2.— Model for the geometry of the HST/NICMOS pupil 
and location of the sample points for the determination of the Ker- 
phase relations. The 156 sample points of the pupil fall on a regular 
square grid with a step of l/16 th of a pupil diameter that do not 
intersect with the central obstruction or the spider arms. 

While not as immediately tangible as the notion of 
closure-phase, this proposed generalization, hereafter ref- 
ered to as kernel-phase (or Ker-phase) since it relates to 
the kernel of the matrix A T , exhibits a unique advan- 
tage over the classical closure-phase: it is not restricted 
to non-redundant apertures and makes it possible to ex- 
tract phase-residual immune information from images ac- 
quired from telescopes of arbitrary pupil geometry. 

This has some obvious advantages over the restrictive 
non-redundant case: 

• throughput: non- redundant aperture 
masks exhibit a typ ical 5-15% throughput 
(Mar tinache et ail l2007t ). and photon noise of 
the companion one tries to detect may be the 
dominant source of noise. Given that it benefits 
from the same phase-noise cancelling properties 
as closure-phase, for a given readout noise and 
exposure time, kernel-phase on an unmasked 
aperture offers an immediate boost in sensitivity 
(or dynamic range) on faint sources. 

• number of observable quantities: a common non- 
redundant aperture mask design exhibits nine sub- 
apertures, therefore forming (|) = 36 baselines 
and ( §) = 28 independent closure phases (|Monnierl 
2000). More independent kernel-phases can be 
extracted from the Fourier transform of a full- 
aperture image, which will provide a better char- 
acterization of the target. 

Another incidental advantage is that, being a product 
of the SVD, all the kernel-phase relations contained in 
K form an orthonormal basis, and therefore do not in- 
troduce correlation in the data. A consequence is that 
manipulating Ker-phases does not require to keep track 
of the covariance matrix used for closure-phases in mask- 
ing interferometry, which simplifies their interpretation. 



2 The use of this model for wavefront sensing purposes will be 
the object of another paper. 
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In discretizing the pupil into a finite number of sub- 
apertures, one important assumption is made: the phase 
(or more generally, the complex amplitude of the electric 
field) is assumed to be uniform within each sub-aperture. 
Yet even for a space borne telescope, in the absence of at- 
mosphere, this is only an idealization as small scale struc- 
tures like polishing imperfections of the primary mirror 
for instance, will impact, to some extent, the value of 
the Ker-phases. This issue is not proper to Ker-phases 
and also affects closure-phases. Thus, unless perfect (i.e. 
single-mode) spatial filtering is performed within each 
sub-aperture of a non-redundant array, the closure phase 
on a point source is not exactly be zero. 

This effect can somewhat be mitigated by substracting 
from the Ker-phases of a science target, the Ker-phase 
signal measured on a point source observed in identi- 
cal conditions. NRM-in terferometry results reported in 
iMartinache et al. (2009) for instance, make extensive use 
of this kind of calibration: from the ground, this ap- 
proach is very powerful as it makes it possible to cal- 
ibrate other sources of systematic errors like the effect 
of broadband filters which smear out the Fourier plane 
and differential atmospheric refraction. From space, this 
may not be as essential depending on the science goal: 
if the Ker-phases obtained on a binary system are non- 
calibrated, then they will simply contain a systematic 
error term that will limit the achievable contrast. 

3. KERNEL PHASE ANALYSIS OF HST/NICMOS DATA 

While the kernel-phase approach may prove difficult to 
apply to ground based observations until extreme Adap- 
tive Optics become available, it can readily be applied 
to diffraction-limited observations made from space. It 
is tested here on a series of non-coronagraphic narrow 
band images acquired with the Near-Infrared Camera 
and Multiobject Spectrometer (NICMOS) onboard the 
Hubble Space Telescope. Two datasets acquired with 
the NICMOS1 in the F190N filter on two distinct ob- 
jects are used: the first target is a calibration star, SAO 
179809, which was observed in 1998; the second is is the 
high-proper motion star GJ 164, around which a compan- 
ion was astrometrically discovered and whose existence 
was confirmed after PSF mo delization and substr action 
of these NICMOS 1 images bv lPravdo etall (|2004f >. This 
latter target is an ideal benchmark: given its expected 
< 10 : 1 luminosity contrast, one should expect a strong, 
unambiguous Ker-phase signal. 

Moreover, ground-based infrared aperture mask- 
ing interferome t ry measurements reported by 
IMartinache et all (|2009f) combined with the astrometry 
have lead to strong constraints on the orbit of the 
companion around the primary. The location of GJ164 
B measured from the Ker-phase analysis of the data can 
be compared to the orbit prediction. 

Figure [2] shows the model of the pupil used for this 
exercise. The HST pupil exhibits a 30 % central obscu- 
ration as well as 90° spider arms (actual dimensions were 
taken from the NIC1 configuration file in the TinyTim 
PSF simulation package for HST). The phase across the 
pupil is discretized into a 156 elementary sub-aperture 
array, whose locations fall on a regular square grid of 
step of l/16 th of the outer pupil diameter. The phase 
sample is assembled into a 155-component (N — 1) vec- 
tor ip. 



These 156 pupil phase samples map in the (u, w)-plane 
onto a square grid of 366 distinct elements 3 . The re- 
sulting (it, w)-sampling is illustrated in Fig. [3] For this 
analysis, A (cf. Eq. [IJ is therefore a 155 x 366 rectan- 
gular matrix, whose SVD reveals that 78 singular values 
are non-zero, leaving 366 — 78 = 288 Ker-phase relations. 

The GJ164 data consists of a total of 80 frames, ac- 
quired at average Julian Date 2453049.3 (February 14, 
2004 UT). Each image is a non-saturated 32 second ex- 
posure, and the target was acquired in a total of 10 dis- 
tinct dither positions. Note that this dataset does not 
include images on a point-source and therefore, the Ker- 
phases calculated from this dataset are non-calibrated. 
Images corresponding to one dither position were simply 
coadded forming a final total of ten 250-second exposure 
images, and assembled into a datacube. The SAO 179809 
data consists of four distinct 20-second exposure frames 
assembled into a separate datacube. 

For both datacubes, the images were then centered 
with sub-pixel accuracy and win dowed by a super - 
Gaussian function as described by iKraus et al.l (|2008l ) 
to limit sensitivity to readout noise. The window size 
is about 25X/D in diameter, which is significantly larger 
than the field of view in which this technique is relevant. 

After this preparatory stage, the images are simply 
Fourier-transformed (cf. second panel of Fig. |3]), and 
the signal R • $ is directly measured for each of the 366 
(u, v) points by sampling the imaginary part of the lo- 
cal complex visibility. The uncertainty associated to the 
measurement of each phase is estimated from the disper- 
sion of the signal in the direct neighborhood of the (u, v) 
point. 

The (u, v) signals are then assembled into Ker-phases 
using the relations gathered in the rows of K and uncer- 
tainties are propagated. The procedure is repeated for 
each of the frames within each datacube. The final re- 
tained series of 288 Ker-phases is the weighted average 
for all frames. 

Because the Ker-phase relations are designed to pro- 
duce quantities independent from pupil phase errors, a 
point source is expected to exhibit zero signal within 
uncertainty. Despite the small number statistics (four 
frames acquired on SAO 179809), the Ker-phase of the 
calibrator do average to zero (with a 19.7° standard devi- 
ation) , while the binary exhibits a large signal amplitude 
(> 100°) in comparison with the uncertainty of individ- 
ual Ker-phases (~ 2°). The third panel of Fig. [3] com- 
pares the Ker-phase histograms of both datasets. 

To further investigate the GJ 164 data, a parametric 
model of the (u, w)-plane phase $o for a binary star is 
needed. The parameters are: the angular separation, the 
position angle of the secondary relative to the primary 
and the luminosity contrast ratio. The model phase $o 
is then multiplied by the diagonal matrix R, and finally, 
transformed into model Ker-phases using the relations 
established during the SVD. 

The agreement between the data and the model is very 
good (cf. panel 4 of Fig. [3]), considering the large number 
of measurements (288) adjuted by only three parameters. 
The uncertainties on the Ker-phases, determined from 
the scatter of the data however lead to a best fit reduced 

3 Note for reference, that a non-redundant array of 156 sub- 
apertures would produce exactly 12,090 distinct (u, v) points. 
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Fig. 3. — From left to right: example of a narrow band (F190N) NICMOS1 image used for this work, visualized with a non-linear 
brightness scale; The Fourier transform of this image. The 366 sample points for the phase in the Fourier domain are overlaid; A comparison 
of histograms of the 288 Ker-phases calculated using the relations identified in Section 12.11 By design, the Ker-phases calculated from 
images of a single star are expected to be zero within uncertainty: the corresponding histogram (gray curve) confirms this expectation. In 
comparison, the Ker-phase histogram of the binary (dark curve) appears significantly larger. The same GJ164 Ker-phases plotted against 
the model of a binary star that best fits the data achieve to convince of the presence of a companion in the data. 

Projected probability density function 

TABLE 1 

KER-PHASE DETECTION OF GJ164B 
IN NICMOS DATA COMPARED TO 
PREDICTION FROM ORBITAL 
PARAMETERS 



Projected probability density function 



Projected probability density function 



Position Angle (deg) 



Fig. 4. — To determine confidence intervals for the parameters 
of the binary, a likelihood analysis comparable to the one presented 
by for closure-phase was performed. These panels show the three 
projections of the likelihood function in the region of best fit. Ex- 
cept for the expected correlation between angular separation and 
contrast ratio for a detection within 1 X/D, the solution is unam- 
biguous and well constrained, demonstrating the elegance of the 
Ker-phase approach. 

X 2 larger than one. A global error term (10°) is then 
added in quadrature to the uncertainty to account for 
a systematic error in the non-calibrated Ker-phase and 
produce a final reduced % 2 = 1. 

One can then proceed with determining the uncer- 
tainty on the parameters of the model fit, by close exam- 
ination of the l ikelihood function, very much like what 
is described in (jMartinache et alj|2009f) . The three pan- 
els of Fig. [4] show the evolution of this function in the 
parameter space region near the best solution. Just like 
with closure-phase data, at angular separations less than 
1 X/D, contrast and separation appear to be correlated. 

The uncertainty on each parameter of the model-fit is 
determined after marginalization of the likelihood func- 
tion over the other two parameters. Despite the noted 
correlation, the constraint on the parameters appears 



Parameter 


Ker-phase fit 


Prediction 


Sep. (mas) 


88.5 ± 3.6 


88.2 


P.A. (degrees) 


100.6 ± 0.3 


100.4 


Contrast 


9.1 ± 1.2 





satisfactory, and the best fit (cf. Table [T]) lies well 
within one a of the position predicted from the orbital 
parameters determined from NRM interferometry from 
the ground. It also matches the location reported by 
iPravdo et al.l (|2004f ). after substraction of a simulated 
PSF from the same data, only with a constraint on the 
position angle improved by a f actor of 10. 

From its (H-K) color index, iMartinache et al.1 (|2009l ) 
were able to conclude that GJ 164 B is of spectral type 
later than M8.5, while the primary is well characterized 
as a M4.5 dwarf. One of the most prominent spectral fea- 
tures for M dwarfs is the broad absorption b and of water 
at 1. 8 jim, getting deeper with later types (pones et al.l 
11994 iLeggett et al.|[200l . The -5:1 contrast ratio 
quoted in the NRM paper was determined over the full 
K s filter (bandwidth 2.0-2.3 p, m). A careful exam ination 
of the spectral sequence by (pones et al.l |1994[ ) reveals 
that for this combination of spectral types, the luminos- 
ity of GJ 164B relative to GJ164A seen in the NICMOS 
F190N filter is expected to drop by 30 to 40 % due to 
the water absorbtion band. The 9 : 1 contrast deter- 
mined from the Ker-phase model (cf Table [1]) in this 
narrow filter reflects this evolution. The analysis of this 
GJ164 data demonstrates the validity of the Ker-phase 
approach, by positively detecting a companion whose ex- 
istence was known beforehand. This < 10 : 1 contrast 
detection was however expected to be easy, despite the 
small angular separation of the detection (0.6 X/D). 

Typical NICMOS1 datasets on a given target usually 
consist of four frames only. The SAO 179809 dataset is 
then a representative example and the statistics of its 
Ker-phase (a — 19.7°) can be used in a Monte-Carlo 
simulation to determine contrast detection limits. 

Because the sampling of the (u, i>)-plane exhibits no 
gap, the sensitivity does not depend on the the position 
angle relative to the central star. One can however ex- 
pect it to be a function of angular separation. A total of 
10,000 simulations were performed per point in the angu- 
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NICMOS data contrast detection limits 
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Fig. 5. — Level of confidence in the detection of a companion 
from the analysis of the HST/NICMOS data with the Ker- phase al- 
gorithm. A darker color indicates a region of lower confidence level. 
Three levels arc highlighted: the 90%, 99% and 99.9% confidence 
levels. At angular separation 0.5 X/D (i.e. 80 mas at A = 1.9/^m), 
a contrast limit better than 50:1 is possible at the 99% confidence 
level. 

lar separation/contrast plane to produce the sensitivity 
map displayed in Fig. [SJ The map highlights the 90, 99 
and 99.9 % confidence level detection thresholds. 

The technique looks promising: for such a dataset, at 
0.5 X/D, a 50:1 contrast detection appears possible at 
the 99 % confidence level. The sensitivity increases and 
peaks at 180 mas, which unsurprisingly corresponds to 
the location of the first zero of the diffraction for the cen- 
trally obstructed telescope (about 1.1 X/D), and reaches 
~200:1. 

4. CONCLUSION 

Classical closure-phase appears to be one special case 
of a wider family of observable quantities that are im- 
mune to phase noise and non-common path errors. In 
the high Strehl regime, it was demonstrated that closure- 
phase like quantities, called Ker-phases, can be extracted 
from focal plane images, and provide high quality "inter- 
ferometric grade" information on a source, even when 



the pupil is redundant. The Ker-phase technique was 
successfully applied to a series of archive NICMOS im- 
ages, clearly detecting a 10:1 contrast companion at a 
separation of 0.5X/D. Non-calibrated Ker-phase appears 
sensitive to the presence of 200:1 contrast companion at 
angular separation 1X/D. Re-analysis of other compa- 
rable NICMOS datasets with this technique might very 
well lead to the detection of previously undetected ob- 
jects in the direct neighborhood of nearby stars. 

Unlike closure-phases, which are extremely robust to 
large wavefront errors, the use of Ker-phases is how- 
ever for now restricted to the high-Strehl regime, and 
will only become relevant to ground based observations, 
when extreme AO systems become available. There is 
nevertheless hope to be able to extend the application 
of Ker-phases to not-so-well corrected AO images, using 
additional differential techniques. One possibility, con- 
sists in using integral field spectroscopy, to follow in the 
Fourier plane, the evolution of the complex visibilities as 
a function of wavelength. With enough resolution and 
spectral coverage, this indeed allows to identify the pha- 
sors contributing to the power contained at one spatial 
frequency. 

The singular value decomposition of the transfer ma- 
trix used to create Ker-phase relations can also be used 
to produce a pseudo inverse to the matrix, and in some 
cases, allows to inverse the relation linking the (u, v) 
phases to the pupil phases. This means that under cer- 
tain conditions, a single monochromatic focal plane im- 
age can also be useful for wavefront sensing purposes. 
This is particularly interesting since the measurement is 
happening at the level of the final science detector, which 
therefore allows to calibrate non-common path errors. 
The application of the formalism to wavefront sensing 
will be the object of a future publication. 



The author thanks Michael J. Ireland and Olivier 
Guyon for the useful discussions of the ideas presented 
in this work. 
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